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Abstract 

The radiative corrections to the pure leptonic decay Ds — ^ii'i up-to one-loop 
, order is presented. We find the virtual photon loop corrections to Dg — ^rvr is 

Qj ' negative and the corresponding branching ratio is larger than 3.51 x 10~^. Consid- 

• ering the possible experimental resolutions, our prediction of the radiative decay 

Dg — >tVt-^ is not so large as others, and the best channel to determine the Vcs 
\ or fo^ is Dg — >^Ufj_'y. How to cancel the infrared divergences appearing in the 

■ loop calculations, and the radiative decay Ds — >£f£7 is shown precisely. It is 



emphasized that the radiative decay may be separated properly and may com- 
pare with measurements directly as long as the theoretical 'softness' of the photon 
corresponds to the experimental resolutions. 

PACS numbers:13.20.He, 12.38.Lg, 12.39.Jh, 13.40.Ks 



1 Introduction 

The pure leptonic decays of the heavy-light meson Dg are very interesting. In principle, 
the pure-leptonic decay Ds — >£f^(see, Fig.l) can be used to determine the decay constant Jd^ 
if the fundamental Cabibbo-Kobayashi-Maskawa matrix element Vcs of Standard Model (SM) 
is known. Conversely if we know the value of decay constant f^s from other method [Q, ^, |3|, 
these process also can be used to extract the matrix element Vcs- But there are the well known 
effect of helicity suppression and there are only one detected finial state, the measurement of 
pure leptonic decays of such pseudoscalar meson is very difficult. 

The pure leptonic weak decays of the pseudoscalar meson corresponding to Fig.l are helicity 
suppressed by factor of / m|)^ : 



Of them only the process Dg — ^ri^r is special, it does not suffer so much from the helicity 
suppression, and its branching ratio may reach to 4.5%. However the produced r will decay 
promptly and one more neutrino is generated in the cascade decay at least, thus it makes 
the decay channel difficult to be observed. For the channels when the lepton is e or and 
their neutrinos, besides the small branching ratios, there are only one detected finial state, the 
measurement of such channels are very difficult. 

Fortunately, having an extra real photon emitted in the leptonic decays, the radiative pure 
leptonic decays can escape from the suppression [Q, ^, ^, furthermore, as pointed out in Ref.0, 
with the extra photon to identify the decaying pseudoscalar meson Dg in experiment from the 
backgrounds has advantages, since one more particle can be detected in the detector. Although 
the radiative corrections are suppressed by an extra electromagnetic coupling constant a, it 
will not be suppressed by the helicity suppression. Therefore, the radiative decay may be 
comparable, even larger than the corresponding pure leptonic decays[§, ^, ^. So, the problem 
to increase the accuracy of the theoretical calculation, at least up to the first order radiative 
corrections, emerges. 

The radiative pure leptonic decays, theoretically, have infrared divergences and will be 
canceled with those from loop corrections of the pure leptonic decays. In all the existing calcu- 
lation of radiative decays ^, |6|, this part is ignored, since they do not include the radiative 
corrections of the pure leptonic decays. In this paper, we are interested in considering the ra- 
diative decays and the pure leptonic decays with one-loop radiative corrections together. Since 
the process Ds — >Ti'r dose not suffer the helicity supression and has a large branching ratio, 
the corresponding loop correction(virtual photon) to these process should has a considerable 
larger branching ratio, at least comparing with the radiative decay, and can not be ignored. 

The paper is organized as follow: in Section II, we present the calculations of radiative 
pure leptonic contributions: i) besides the leptonic decays shown as Fig.l, the radiative decays 
into three family leptons with a real photon appearing in the final state as shown in Fig. 2; ii) 
the virtual photon corrections to the pure leptonic decays, i.e., a photon in loop as shown in 
Figs. 3.1, 3.2, 3.3, 3.4. In Section III, we evaluate the values of the decays. As uncertainties the 
dependence of the results on the parameters appearing in the considered model is discussed. 
Finally, preliminary conclusion is obtained. 



Let us start with the radiative decay. These contributions are corresponding to the four 
diagrams in Fig. 2. According to the constituent quark model which is formulated by Bethe- 
Salpeter (B.-S.) equation, the amplitude turns out to be the four terms Mi{i = 1,2,3,4): 



2 The Radiative Leptonic Contributions 
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where <?) is Bethe-Salpeter wave function of the meson Dg; p is the momentum of D^; e, 
are the polarization vector and momentum of the emitted photon. In the quark model, the 
momenta of s, c-quarks inside the bound state, i.e., the Dg meson, are: 
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where q is the relative momentum of the two quarks inside the Dg meson. As the Dg meson 
is a nonrelativistic bound state in nature, so the higher order relativistic corrections may be 
computed precisely, but being an approximation for a iS-wave state, and focusing the light on 
the radiative decay corrections only at this moment now, we ignore q dependence, i.e., we may 
still have the non-relativistic spin structure for the wave function of the meson Dg (a state) 
correspondingly: 

Here V'(O) is the wave function at origin, and by definitions it connects to the decay constant 
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where m is the mass of Dg meson. Moreover we note that for convenience we take unitary 
gauge for weak bosons to do the calculations throughout the paper. With a straightforward 
computation, the amplitude can be simplified as: 
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where ki is the momentum of the charged lepton, and 
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As a matter of fact, there is infrared infinity when performing phase space integral about 
the square of matrix element at the soft photon limit. It is known that the infrared infinity 
can be cancelled completely by that of the loop corrections to the corresponding pure leptonic 
decay Ds tv. 

In Eqs.(7) and (8), the infrared terms can be read out: 
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As the diagrams (g), (h), (i), (j) in Figs. 3. 3, 3.4 always have a further suppression factor 
rr? jm^ to compare with the other loop diagrams, and there is no infrared infinity in these four 
loop diagrams, we may ignore the contributions from these four diagrams safely. Furthermore 
we should note that in our calculations throughout the paper, the dimensional regularization 
to regularize both infrared and ultraviolet divergences is adopted, while the on-mass-shell 
renormalization for the ultraviolet divergence is used. 

If Feynman gauge for photon is taken (we always do so in the paper), the amplitude 
corresponding to the diagrams (a), (b) of Fig.3.1 can be written as: 
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where the I, k2 denote the momenta of the loop and the neutrino respectively. These two terms 
also have infrared infinity when integrating out the loop momentum I. 

After doing the on-mass-shell subtraction, the terms corresponding to vertex and self-energy 
diagrams (c), (d), (e), (f) can be written as: 
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Now let us see the cancellation of the infrared divergencies precisely. The infrared parts 
of the decay widths which are from the interference of the self-energy and vertex correction 
diagrams with the tree diagrams: 
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(13) 



the infrared part of the decay widths from the interference of the "box" correction diagrams 
with the tree diagrams: 
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the infrared part of the decay width from the real photon emission: 
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Here A£^ is a small energy, which corresponds to the experimental resolution of a soft photon 
so that the phase space of the emitting photon in fact is divided into a soft and a hard part by 
IS.E. Where /U is the dimensional parameter appearing in the dimensional regularization. It is 



easy to check that when adding up all the parts: the real photon emission 5T 



virtual photon corrections ^p^W'^^red^ ^p«^/''afed^ ^j^^ infrared divergences ^ — 7 are canceled 
totally and the /i dependence is also cancelled. Hence we may be sure that we finally obtain 
the pure leptonic widths for the Dg meson decays to the three families of leptons which are 
accurate up-to the 'next'-order corrections and 'infrared finite' but depend on the experimental 
resolution IS.E. 
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3 Numerical Results and Discussion 

First of all, let us show the 'whole' leptonic decay widths, i.e., the sum of the corresponding 
radiative decay widths and the corresponding pure leptonic decay widths with radiative correc- 
tions, and put them in Table (1). Why we put the radiative decay and the pure leptonic decay 
with radiative corrections together here is to make the width not to depend on the experimen- 
tal resolution for a soft photon. The values for the parameters a = 1/132, \Vcg\ = 0.974[^], 
mn, = 1.9686 GeV, mg = 1.5 GeV, rric = 1.7 GeV and fo, = 0.24 GeV[|. 



Table (1) The 'Whole' Leptonic Decay Widths (in unit GeV) 



Te{W-") 


3.593 


r^(io-^5) 


6.604 


r,(io-i4) 


5.807 



For comparison, the width of each pure leptonic decay at tree level with the same param- 
eters as those in Table (1) is put in Table (2). 

Table (2) The Pure Leptonic Decay Widths (in unit GeV) of Tree Level 



re(io-^^) 


1.521 


r^(i0"^5) 


6.463 


r,(io-i4) 


6.30 



If the lifetime of Ds meson is t{Ds) = 0.469 x 10 the corresponding branching ratios 

are showed in Tables (3), (4). 

Table (3) Branching Ratios of the 'Whole' Leptonic Decays 



Se(10-^) 


2.56 




4.706 


5r(10-^) 


4.138 



Table (4) Tree Level Branching Ratios of The Pure Leptonic Decays 



Se(10-0 


1.084 




4.605 


i?r(10-^) 


4.489 



We can see that, the whole decay widths Tg and T^ in Table (1) are larger than the 
corresponding decay widths of tree lever in Table (2), while the Tr in Table (1) is smaller than 
the one in Table (2). This means the sums of radiative leptonic decay and loop correction 
as the first order contributions of pure leptonic decay are positive for the decays with a e or 
/i in the finial state, and negative for the one with a r in the finial state. It also show that 
the contributions of loop corrections are negative, the dominate contributions of first order 
corrections to the pure leptonic Ds decays are radiative decays when the lepton is e or /i, 
and is loop corrections when the lepton is r. So, the loop contributions are important for the 
decays Dg — > //i/^ and Dg — > tv^, especially for the later. Through Table (3) and (4), we 
obtain that the rediative decay has a branching ratio Br{Ds IJ-^til) > 1-01 x 10"'^ and the 
loop correction to Dg tv-j- has a branching ratio Br > 3.51 x 10^"^. 

To see the contributions of the radiative decays precisely we present the radiative decay 
widths with a cut of the photon energy, i.e., the widths of the radiative decays Dg — > Ivy 
with the photon energy E.y > kmin as the follows: kmin = 0.00001 GeV, kmin = 0.0001 GeV, 
kmin = 0.001 GeV, kmin = 0.01 GeV and kmin = 0.1 GeV respectively in Table (5). 

Table (5): The Radiative Decay Widths (in unit GeV) with cuts of the photon momentum 
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0.00001 


3.582 


6.878 


8.892 


0.0001 


3.581 


5.484 


6.452 


0.001 


3.580 


4.091 


4.025 


0.01 


3.578 


2.704 


1.708 


0.1 


3.474 


1.363 


1.020 



The corresponding branching ratios of the radiative decay with cuts of photon momentum are 
showed in Table (6) 

Table (6): The Radiative Decay branching ratios with cuts of the photon momentum and the 

results of Ref g |, | 
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Bre 
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GeV 


10"^ 


10-4 


10-** 


0.00001 


2.552 


4.901 


6.336 


0.0001 


2.552 


3.908 


4.597 


0.001 


2.551 


2.915 


2.868 


0.01 


2.549 


1.927 


1.217 


0.1 


2.475 


0.971 


0.727 


Refi 


10 


1 




Refi 


17 


1.7 




Refi 


7.7 


2.6 


320 



Considering the possible experimental resolutions, our prediction of the radiative decay 
widths T{Ds — > evej) and T{Ds — > fJ-i^^^) are close to the values in RefQ, ^, ^, but our 
prediction of T{Ds tVtI) is much smaller than the one in Ref|^]. In our model, if we using a 
smaller cut kmim the obtained Br[Ds — > Ivi'^) will become larger, but since the decay widths 
depend on Log{kmin), the change of branching ratios will be not so much on the selecting of 
kminj we can see this in Table (5) and (6), and for another example, if k^nin — 

1 X 10"^°, the 

Br{Ds evel) = 3.59 x 10-^ Br{Ds nu^j) = 1.38 x lO^^, Br{Ds rvrl) = 2.11 x 10-^ 
but it is very difficult to obtain so small a k^in in experiment. We can conclude that the best 
radiative decay channel is easy to search is Dg — > ^i^^'f- 

For the convenience to compare with experiments, we present the photon spectrum of the 
radiative decays in Fig.4 and Fig.5. In addition, we should note that the widths are quite 
sensitive to the decay constant , and are sensitive to the values of the quark masses and 
rric- 
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Figure 3: 1. Box- loop diagrams for Dg 




Figure 3: 4. Self-energy diagrams for Ds — > ii^. 
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ure 4: Photon energy spectra of radiative decays Dg — > ii^e'yii — 
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Figure 5: Photon energy spectra of radiative decays D, 



